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The fusion of the probabilistic finite element method (PFEM) and reliability analysis 
for probabilistic fracture mechanics (PFM) is presented. The authors have developed a 
comprehensive method for determining the probability of fatigue failure for curved crack 
growth. The criterion for failure or performance function is stated as: the fatigue life of a 
component must exceed the service life of the component; otherwise failure will occur. An 
enriched element that has the near crack-tip singular strain field embedded in the element is 
used to formulate the equilibrium equation and solve for the stress intensity factors at the 
crack-tip. The crack growth law (e.g., Paris law) and crack direction law (e.g., direction 
which gives rise to the maximum tangential stress at the crack-tip) are also incorporated into 
the formulation. The loading is mixed-mode (i.e., opening mode and in-plane shear mode 
which are designated mode I and mode II, respectively) with randomness in the applied 
loads, material properties, component geometry, crack geometry, crack growth law, and 
the crack direction law. Examples of random crack geometry are the initial and final crack 
lengths, initial crack angle, and initial crack position. The methodology consists of 
calculating the reliability index, which is used to calculate the first-order probability of 
failure, by solving a constrained optimization problem. Included in the optimization for 
reliability is the equilibrium equation for stress intensity factors, the crack growth law, and 
the crack direction law. 

Performance and accuracy of the method is demonstrated on a classical mode I 
fatigue problem. The initial and final crack lengths, material parameters in the Paris law, 
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and the applied stress are modeled with random variables. The results imply that the 
service life of a component must be well below the deterministic fatigue life of the 
component when many parameters are considered random. The fusion of PFEM and 
reliability for fatigue crack growth is computationally quite efficient and provides a 
powerful tool for the design engineer. 
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SCOPE OF WORK 


Probabilistic Finite Element Method fPFEMl 

• Uncertainties 

Random Load 

Random Material Properties 
Random Geometry' 

• Statics and Dynamics 

• Elasticity and Nonlinear Elasticity 


Brittle Fracture bv PFEM 

• Probability of the Stress Intensity Factor exceeding 
the Fracture Toughness of the Material 


P ( K i > K k ) 

P ( G > G c ) 


Fatigue Crack Growth bv PFEM 

• Probability of the Desired Service Life exceeding the 
Fatigue Life 


P ( T s > T ) 


Introduction 


Probabilistic Fracture Mechanics fPFM ) 

• Model the crack-tip singularity in stress. 

• Determine statistics of stress intensity factors. 

• Probability of fracture and fatigue failure. 


Uncertainties which need to be considered in PFM are: 

• Random load, material properties, and geometry; 

• Random crack geometry; 

length, 

location, 

orientation, 

• Random crack growth laws; 

• Random crack direction laws; 

• Random correlation among detectable cracks and 
micro-cracks. 
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Reliability of Components 
Subject to 

Fatigue Crack Growth 
by 

Probabilistic Finite Elements 


Methodology and Objective 

• Given : 


• Component with Crack Subject to Cyclic Stress 

• Method for Determining Stress Intensity Factors 

• Crack Growth Law 

• Crack Direction Law 

• With Randomness in : 

• Applied Loads 

• Component Geometry 

• Crack Geometry 

• Crack Direction Parameters 

• Crack Growth Parameters 

• Material Properties 

9 Determine Probability of Fatigue Failure 

• Probability that the desired service life is greater 

than the fatigue life. 

P ( T s > T ) 


Stress Intensity Factor 



Kj as mode I stress intensity factor 
Kjj = mode II stress intensity factor 
K III = m0( ^ e m stress intensity factor 


Consider only Modes I and II: 
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Random Variables 


System Uncertainties 

• Applied Loads 

• Component Geometry 

length, height, width, etc. 

• Crack Geometry 

length, orientation, position 

• Crack Direction Parameters 

• Crack Growth Parameters 

• Material Properties 

elastic modulus, poisson’s ratio, 
fracture toughness, etc., 


Modeled by a correlated random variable vector 
b (q x 1) 


Transform b to a Independent Standard Normal Random 
Variable Vector 

r (qx 1) 


* Eigenvalue Orthogonal izati on (Normal R.V.'s) 

Covariance Matrix for b 

• Rosenblatt Transformation (Non-normal R.V.’s) 

Joint Probability Distribution for b 


Reliability Analysis 


Performance Function or Margin of Safety 

g>0 g = 0 g < 0 

safe state limit-state failure state 

Probability of Failure 

p f = P(g ^ 0) 


First-Order Estimate of the Probability of Failure 
P fl = <£ (-P) 

O is the standard normal cumulative distribution, 
P is the reliability index. 


Reliability Index 

Minimum Distance from Limit-State Surface to 
Origin in Independent Standard Normal Space 

p 2 = <^7 
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Crack Growth Law 


Experimental fatigue crack growth data for different 
materials. 



ir. ■ l 

AK e q is equivalent stress intensity factor range 



Region A: A sufficient level of activity at the crack-tip has 
not been attained resulting in very slow crack 
growth. 

Region C: Corresponds to rapid unstable crack growth. 

Region B: Represents the zone where most fatigue crack 
growth is analyzed. 


Paris Law: ^ = D (AK eq ) n 

D is a material parameter 
n is a material parameter 
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Integrating from the initial crack length to the final crack 
length gives in the fatigue life 


T = 


C da 
D (AK eq ) n 


Equivalent Stress Intensity Factor 
* Maximum Stress Theory 


K eq = 0 x 


2\Q 

0 = COS^ 


2|_ 

12 


COS 


3 sin 


CRACK DIRECTION LAW 


• Maximum Tangential Stress Criterion 

• Minimum Strain Energy Density Criterion 

• Maximum Energy Release Rate Criterion 


Maximum Tangential Stress Criterion 

* Crack grows in a radial direction from the crack 

• Crack grows in the plane perpendicular to the 
direction of maximum tension 


Crack grows in the direction where the tangential 
stress is a maximum and the shear stress is zero. 

Z(tc, 0) = <J) T x = 0 

f sin 0 1 


= 


L 3 cos 9 - 1 J 

For mode I fracture, 9 = 0° and is constant 
For mode II is fracture, 9 = 70.54° and changes 
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PERFORMANCE FUNCTION FOR FATIGUE 


Basis for determining the reliability and/or probability of 
failure depends on the choice of a performance function. 


P f = P(g £ 0) 


Performance function for fatigue crack growth 


T $ is the desired service life 
T is the fatigue life 

The component fails when the fatigue life is less than 
the desired service life. 


DISCRETIZATION 


Discretize the Crack Path into "nint" Integration Points 

• Crack Length 

a i = 1 t(a f - aj)£j + (a f + a,)] i = 1 nint 

• Crack Growth Law 


T = 


+ 1 

r J 


- 1 


D (K eq ) n 


• Crack Direction Law 


J = 


da 


^k = ^k x k = ® k = 1, .... nint 

Crack direction is calculated at each integration point 
Crack grows to the next integration point 

• Equilibrium 
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RELIABILITY AS AN OPTIMIZATION PROBLEM 
FOR FATIGUE CRACK GROWTH 

First-Order Probability of Failure 
-4 Reliability Index 

-4 Constrained Optimization Problem 

Minimize the Reliability Index 



Subject to Equality Constraints 
Crack growth direction law 

= 0 k = 1, nint 

Equilibrium 

5- = F- i = 1, nint 

Subject to Inequality Constraint 

Limit-State Surface and Failure State 
T - T s < 0 
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Define a Lagrange Functional 


& ( b > ^1 ^nint- ®1» 5 ninf a - 

•Pi-- (Point- 9 1 e nint> 

nint 

= p 2 + X m- T [ F i - K i 6 i 3 

i=l 

nint 

+ X *Pk £ Z k = 0 1 

k=l 

+ ^ {t • T s + a 2 } 


Necessary Conditions 

(derivatives w.r.t. independent variables) 

Probabilistic Finite Element Method 
(PFEM) 

Eliminate Lagrange Multipliers 


Final Form of Optimization for Reliability 


+ X L g = 0 


L P 


ap 2 


“ db 


nint I 

r i 

' T ’ 9 k 1 


9K: 

ii 

X 

k=l 

C Zk )e k (Zk) ", 

+ T ’K- 

1 

db 

L . \ 

nint 

T -e. i 



+ 

X 

k=l 

We k (2kH ; 

+ T b 



Solve for b via an Optimization Algorithm 

* Lagrange Multipliers 

* Gradient Projection Algorithm 

• HL-RF Algorithm (Tangent Matching) 

• Modified HL-RF Algorithm 
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CRACK GROWTH MOVING MESH 


• Variation of ALE formulation 

Pseudo-Static (neglect transport term) 

* Moving Element Mesh with Constraints 

Aspect ratio of elements must remain one 


b 
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EXAMPLE 


Single-Edged Cracked Plate with an Applied Load 
b = [ a } , a f , D , n , a ] T 

aj is the initial crack length 
a f is the final crack length 

D is a material parameter 
n is a material parameter 
a is the applied stress 


Parameter 

Mean 

Standard Deviation 

Percent 

Lengrn (Lj 

10.0 in 

0.0 

0.0 

Width (W) 

4.0 in 

0.0 

0.0 

Thickness (t) 

1.0 in 

Y. 

0.0 

Young's Modulus (H) 

30,000.0 ksi 

0.0 

0.0 

Poisson’s Ratio (v) 


0.0 

0.0 

Applied Stress (a) 

12.0 ksi 

3.0 ksi 

25.0 

Initial Crack Length Cap 



100.0 

Final crack Length (ap 

0.1 in 

0.01 in 

10.0 

Fatigue Parameter (D) 


3.0 x 10° 1 

30.0 

Fatigue Parameter fn) 

3.25 

0.08 

2.5 


• Variations in the parameters are reasonable. 

• nint = 50 

• Deterministic fatigue life, T = 7,500,000 cycles 
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Reference Solution 

Kj = 'y a Va 

^ = 1-99 - 0 . 41 ^+ 18 . 7 (^ r ) 2 - 38 . 48(^) 3 + 53 . 85 (| r ) 4 


Deterministic Stress Intensity Factor for Reference and FEM Solutions 
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Deterministic Fatigue Life 
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Probability of Failure for Combined Randomness 
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